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Abstract. In this paper we analyze the non-integrability of the Wilbeforce 
spring-pendulum by means of Morales-Ramis theory in where is enough to 
prove that the Galois group of the variational equation is not virtually abclian. 
We obtain these non-integrability results due to the algebrization of the vari- 
ational equation falls into a Heun differential equation with four singularities 
and then we apply Kovacic's algorithm to determine its non-integrability. 



1. Introduction 

The Wilberforce pendulum consists of a mass hanging on a flexible spiral spring 
that is free to oscillate in both the standard longitudinal mode and the torsional 
mode. When the mass is lifted above its equilibrium point and released from rest, 
it oscillates up and down along a vertical line, slowly transferring its energy into a 
rotational oscillation. If the nuts screwed onto the vanes protruding from the sides 
of the mass are adjusted to give the appropriate moment of inertia such that the 
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frequencies of the longitudinal mode and the torsional mode are the same, the pen- 
dulum will transfer its energy back and forth completely between these two modes 
os oscillation, see [3]. We shall refer to the Wilberforce spring-pendulum when the 
spring is free to swing in a plane, thus adding an extra degree of freedom. This sys- 
tem displays evidence inherent to chaotic dynamical systems. Hence, establishing 
its non-integrability is our goal. The main tool for studying the non-intcgrability 
of this kinds of Hamiltonian systems is Morales-Ramis theory. One important cri- 
terion to obtain non-integrability of autonomous Hamiltonian systems by means 
of differential Galois theory is due to J. J. Morales- Ruiz and J.P Ramis [10], and 
Morales- Ruiz [9]. The Morales-Ramis theorem connects two notions: integrability 
of Hamiltonian systems and integrability of linear systems of equations. In partic- 
ular, Morales-Ramis theory has been used to prove the non-integrability of spring 
pendulum systems, see [5, 7, 9]. 

This paper is organized as follows. We start with a brief description of the 
basics of the Morales-Ramis theory of non-integrability of the Hamiltonian sys- 
tem, geometric objects associated and the theory of linear homogeneous system 
of differential equation with singular points, following [4] and [1]. In the section 
4 we present the necessary theoretical background to apply Morales-Ramis theory 
to the Wilberforce spring-pendulum Hamiltonian system. Finally, the Section 5 
is devoted to a detailed analysis of apply the Theorem 3.4 in order to prove that 
the Wilberforce spring-pendulum Hamiltonian system is non-integrable in terms of 
meromorphic first integrals. 

2. Variational equations and connections 

Let M be complex analytic manifold and X a meromorphic vector field on M. 
Let us denote by K the field of meromorphic functions on M and d: K — >■ K the 
derivative with respect to X so that (K, d) is a differential field. From now on we 
will write /' for df with / an element of K. 

Let 7r: E — > M be a vector bundle, that we assume to be meromorphically 
trivial. Let us recall that this always holds in the algebraic case due to the vanishing 
of the first Galois cohomology group of the linear group, and also holds if M is an 
open Ricmann surface, due to Brauer theorem. 

Let E be the space of meromorphic sections of tt. A linear connection in the 
direction of X is a K d linear map V : E — >• E satisfying Leibniz rule V(/e) = 
(df)e + /V(e). Let us consider two vector bundles (E, V) and (E', V') over M 
endowed with linear connections in the direction of X. A morphism of connections 
is a AT-linear map (j>: E — >• E' such that <j) o V = V o (f>. 

Let {ei, . . . , e„} be a basis of E. The n 2 coefficients a,j £ K in the following 
expression arc uniquely determined. 

n 

V(e t -) = y^ajj-ej, i = l,...,n 

3 = 1 

Therefore the so-called connection matrix A = (a^) £ Mat(n x n,K) determines 
uniquely V. 

A meromorphic section e £ E is called horizontal if it satisfies V(e) = 0. Let 
us write e = Y17=i x i e ii x = i x i> ■ ■ ■ > x nY an d x' = 8x whenever c <S K s . It follows 
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from equation (1) that x satisfy the system of linear differential equations, 

x' = -A'x. 

Therefore, the choice of a linear basis of E over K gives a one-to-one correspondence 
between linear connections in E and n x n systems of linear differential equations 
with coefficients in K. 

Remark 1. Let us assume now that M is a Riemann surface T, the differential 
Galois theory in this case is developed in [15], and can be generalized to the case 
in which the held C of constants of K is the field of complex numbers C, id est, the 
vector field A does not admit meromorphic first integrals. To a connection V in the 
direction of A it corresponds an algebraic group Gal(E, V) which is embedded into 
GL(n, C) up to a conjugacy class. In the development of our application we will 
just make use of the following geometric properties of the differential Galois group. 

(1) An exhaustive morphism of connections ir: (E, V) — > (E, V) induces an 
exhaustive group morphism ir* : Gal(E, V) — > Gal(E, V). 

(2) Let us consider (E, V) = (Ei E2, Vi © V2). Then the canonical projec- 
tions induce an injective morphism of algebraic groups, 

7Ti* x 7r 2 „ : Gal(E, V) -> Gal(Ei, V x ) x Gal(E 2 , V 2 ). 

Those geometric properties and many other were exhaustively studied in the more 
general case of the Lie-Vessiot systems in [4]. The necessary differential Galois 
theory background here is given in Appendix A. 

2.1. Generic variational connection. Let us consider M, X, K and d as 

in Section 1. There are two equivalent ways of introducing the variational equation 
for X. The vector field X can be naturally prolonged to the tangent bundle TM 
by application of the chain rule as done in [1]. Equivalently we can consider the 
Lie bracket with A as a linear connection. Let TM be the space of meromorphic 
vector fields in M, then 

: TM -> TM, Y 1 — y V^F = [A, Y] 

is a linear connection in TM in the direction of A. This is called the first variational 
connection of A on the generic point of M 

2.2. Normal variational connection and equations. Let N be an invari- 
ant manifold of M with respect to the flow of A. The following exact sequence of 
vector bundles over N, 

-> TN -> TM| N -> Norm(N, M) -> 

implicitly defines the normal bundle to N in M. It is clear that the connection 
restrict to TM|n and TN and therefore we have also an exact sequence of 
connections. 

->■ (TN,V f | N ) -> (TM| N ,V^| N ) -)■ (Norm(N, M), V) -> 

Let us consider a local system of coordinates x±, . . . , x r , y±, . . . , y m in M such that 
the equations of N are written, 



yi = 0, . . . , y m = 0, 
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and those of X, 



i=l 0X1 1=1 " Vl 

Then, the equations of horizontal section for V, the variational equation along X 
restricted to N arc written, 



(1) 



3 = 1 



N 



N 



m a 



N 



In this system of coordinates, we can easily project the equations onto the 
normal bundle obtaining the equations for horizontal sections of V, 



(2) 



H = E 



dv 3 



N 



Let r be an invariant curve (Ricmann surface) of X contained in N. Its invari- 
ance by X means that the connections discussed above specialize to the bundles 
restricted to T, 

(3) -> (TN| r , V^ |N |r) -> (TM| r ,V_ ? | r ) ->■ (Norm(N, M)| r , V|r) 

The connection is called the variational connection of J? along the integral 
curve F, and V|r is called the variational connection of X normal to N in M along 
r. The equations of horizontal sections for this last one are written 



(4) 



a square system of m linear differential equations with meromorphic coefficients in 
r, or cquivalently a meromorphic linear connection in Norm(N, M)|r- This is the 
so called variational equations of the flow of X normal to the invariant manifold N 
along the integral curve T. 

Remark 2. It is important to remark that the invariant curve T may not 
corresponds to a unique trajectory of the flow of X but also contain equilibrium 
points. This is important because equilibrium points in T are poles of (X|rT) _1 
and therefore singularities of the system (4). In many applications we also replace 
the original manifold M by some suitable compactification such that the considered 
invariant curve T has some equilibrium points at the infinity (see [9] pp. 70-74). 
The curve T may then contain singular points like nodes or cusps. This is coherent 
with the computations above, providing that the field of meromorphic functions on 
a singular Ricmann surface coincides with the one of its desingularization 



3. Non-integrability criterium 



3.1. Morales-Ramis theory. From now on let M be a complex symplectic 
manifold, of dimension 2n with simplectic form fl. Let us consider a symplectic 
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system of coordinates x\, . . . , x n , j/i, . . . , y n . If H is an holomorphic function in M 
then the Hamiltonian vector field Xh is defined intrinsically by the formula, 

dH + i 9 il = 0, 

Xh 

and its expression in local coordinates is, 

dxi dH dyi dH 
dt dyi ' dt dxi 
The symplcctic structure induces a Poisson bracket for functions defined on M, 
namely {F, G} = XpG. The functions F, G in M are said to be in involution if 
{F, G} — 0. A Hamiltonian system is called completely integrable by meromorphic 
functions if it admits n independent meromorphic first integrals in involution. Let 

r be an invariant curve for Xu- Let us define the variational connection V-s in 

x H 

direction of X as in Section 1. Then , (TM|r, V^ h |r) is the variational conection 

to X along the integral curve F. The equation of horizontal sections, so called the 
first variational equation along F is written in symplectic coordinates, 

™ / d 2 H dH \ 

(5) 3 

( d 2 H dH \ 

The Morales-Ramis theory, relates the Liouville integrability of Hamiltonian 
system and the Galois groups of differential equations. The following result (Morales- 
Ramis [10]) and subsequent generalizations are some of the most effective and useful 
known for the theoretical proof of non-integrability of Hamiltonian systems. Let 
us recall that an algebraic group G is said to be virtually abelian it its connected 
component of the identity G° is abelian (see Appendix A) . 

Theorem 3.1. Let H be a Hamiltonian, and Xh its associated Hamiltonian 
vector field in M. Let T be an invariant curve for H. Assume that H is completely 
integrable by meromorphic functions independent in a neighborhood of T but not 
necessarily on T itself. Then the group Gal(TM|r, V^ h |r) is virtually abelian. 

3.2. Analysis of normal variational connection. From now on let N be a 

symplectic submanifold of M, i.e., the pair (N, is a symplectic manifold. Let 
2r be the dimension of N and n = r + m. 

Lemma 3.2. There is a natural decomposition of the bundle TM|n ->N os the 
direct sum of the tangent bundle to N and the normal bundle to N in M, 

TM| N = TM © Norm(N, M). 

PROOF. Let us consider TN 1 C TM the bundle of vectors orthogonal to TN 
with respect to the symplectic form Q. By hypothesis £1\n is non-degenerated and 
therefore TN 1 - is a supplementary bundle for TN. The exact sequence (3) identifies 
TN- 1 with NorTO(N,M) and we get the result. □ 

Lemma 3.3. The variational connection V ' ^ H In in TM|n splits as direct sum of 
the variational connection of V if , in TN and the normal variational connection 

J X H In 

V in Norm(N, M) , 

| N = V? , ©V. 

Xh' a h |n 
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Proof. It suffices to proof that TN 1 - is an invariant bundle for |n- Let 
us take Y meromorphic section of TN^ and Z meromorphic vector field in M. Let 
us note that £^£1 = f2 for any Hamiltonian vector field, and that [X, Z] is defined 
as vector field tangent to N. Therefore, 

CL(V j t\ N Y,Z)=n([X,Y],Z) = (C^Q)(Y,Z)-Q(Y,[X,Z]) = 

and V^|n? is orthogonal to TN with respect to SI. □ 

Theorem 3.4. Let H be a Hamiltonian system in M. Assume that N is invari- 
ant by the flow of Xh- Assume that Xh\n is completely integrable, as Hamiltonian 
system in N, by meromorphic functions. Let T be an invariant curve of Xh in 
N. Then the differential Galois group of the variational equation to the flow of Xh 
along T, Gal(V^|r) is virtually abelian if and only if the differential Galois group 
of the variational equation to the flow of Xh normal to N along T, Gal(V| r ), is 
virtually abelian. 

Proof. First, let us assume that Gal(Vjj|r) is virtually abelian. We consider 
the exact sequence (3). By Remark 1 there is exhaustive group morphism, 

7r 2 : Gal(V^|r) -> Gal(V| r ). 

Therefore, if GalCV^Ir) is virtually abelian then its quotient Gal(V|r) is also 
abelian. Second, let us assume that Gal(V|r) is virtually abelian. From Lemma 3.3 
we have V^|r = V^| N |r © V|r, and by Remark 1 (2) there is an injective group 
morphism, 

(tti x 7r 2 ): Gal(V^|r) -> Gal(V^ |N |r) x Gal(V|r). 

Assume that Gal(V|r) is virtually abelian. From Theorem 3.1 we have that 
Gal(Vj^i |r) is virtually abelian. It follows that Gal(Vj^) is virtually abelian. □ 

The following corollary directly follows from Theorem 3.1. Therefore, we can 
ask what is the use of Theorem 3.4. There are examples, e.g. [12], of application of 
the Moralcs-Ramis theorem in which the non-integrability is not seen in the normal 
variational equation but in the total variational equation. Theorem 3.4 describes 
a theoretical situation in which the non-integrability, if captured by the first order 
variational equation then it is captured with the normal variational equation. 

Corollary 1. Let H be a Hamiltonian system in M. Assume that N is invari- 
ant by the flow of X^j. Assume that Xh\n * s completely integrable, as Hamiltonian 
system in N, by meromorphic functions. Then if Xh is completely integrable in 
M by meromorphic functions then the differential Galois group of the variational 
equation to the flow of Xh normal to N along T, Gal(V|r), is virtually abelian. 

4. The Wilbeforce-spring pendulum 

The classical Wilbeforce spring-pendulum consists of a solid cylinder attached 
to a spring. The cylinder is free to perform oscillations along the vertical. At the 
same time when the spring is elongated it produces a torque proportional to the 
angle of torsion of the cylinder 4> (see Figure 1). We add an extra degree of freedom 
by allowing the spring to swing on a fixed plane in space containing the vertical 
motions. Let r, 9 denote polar coordinates, where r is the distance from the fixed 
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point of the spring to center o mass of the cylinder and 9 the angular deviation 
from the vertical axis. The Lagrangian is 

(6) (r 2 + r 2 9 2 ^j + ^| ifc(r - r ) 2 - ^ i(r - £)e<j) + gmrcosO 

where k is the constant of the spring, tq is the unstretched length of the spring, A 
gives the constant of proportionality of the angular restoring and e is the coupling 
among the vertical and torsional modes. The mass of the cylinder is m and J its 
moment of inertia. 

//////////////////////// 




Figure 1. The Wilberforce spring-pendulum. Taken from refer- 
ence [3]. 



Remark 3. Observe that for e = we recover the spring-pendulum system 
whose aspects of non-integrability have been studied among others by [5], [7], [9] 
and [8]. 

Remark 4. The invariant subsystem 9 — 0,9 = corresponds to the classical 
Wilberforce pendulum which is a quadratic Lagrangian thus a linear system. 

4.1. Adimensionalization. It will be convenient to reduce the seven pa- 
rameters m, J, fc,ro,e,5 to four by the introduction to variables a follows: Let 
£ = ?'o + mg/k, lu 2 = g/£, lu 2 = k/m. Observe that 

r mg = r / ujA 2 = ro 
I H I \uj t 1 

where the parameter 

'-&)' 

is the ratio of frequencies of the linear pendulum and the spring. Then if r = £p, 



r-r =£(p-l + f). 
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Also we introduce a new time variable r via 4r = w s , with this substitution the 
Lagrangian (6) becomes 

L = W(p' 2 +p 2 9' 2 ) + ^0' 2 - i^ 2 (p- 1 + /) 2 - ^ 2 - - l)J> + gmlpco S 6 
I Lm I 11 

a further factorization yields 



L = M 2 



1 e . . am 

where prime (') indicates differentiation with respect r. 

One checks easily that the constant factor can be dismissed leading to the 
adimensional Lagrangian, 

I -i^ ,+ « + I5>*'-i<'- 1 ^-i^-iB^"* + ¥"" t 

Let us introduce the following (dimensionless) parameters 

J A e 

then the Lagrangian finally becomes 

(7) L = i( P ' 2 + pV 2 ) + \a$ 2 - \{p - 1 + /) 2 - ife0 2 - ±c(p - 1)0 + /pcostf. 
The Hamiltonian corresponding to this coordinates assume the form 

(8) H = \ (Pp+y) +^ + l(p-l + f) 2 + lb^ + ^c(p-l)^-fpcose. 

In fact by scaling the angle <f> and its conjugate momenta we can suppose in what 
follows that a = 1 resulting with three parameters: b, c, / 

(9) H=\ (Pp+^f+Pl) +\(P- l + f) 2 + \w 2 + \c{ P - 1)0- fpeosO. 
4.2. The Wilberforce linear subsystem. The equations of motion are given 

by 



(10) 



pf 


= P P , 


e' 


Pe 




p 2 


0' 


= P<j>, 




Pe 


p' 




p 


P 3 " 


p' 

1 6 


= ~fp 


P' 

1 <b 


= -b<j> 
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From these equations, one checks easily that 6 = Pg = is an invariant subsys- 
tem, as asserted previously. The systems reduces in this case to the linear system 



Pi, 



(11) 



c 



Pp = -p+l-^-r, 

That can be recast as two coupled oscillators 

(12) p" + p=l- C -^ 

0" + ty = -fp+f(l-/) 
with Hamiltonian (set Pg = 6 = in (9)) 

(13) #o,o = \ {P 2 P + Pi) + \{ P -l + ff + \btf + \c{p - 1)0 - fp. 



Let us exhibit the reduced system (12) as an integrable Hamiltonian system in a 
convenient normal form. For this purpose let z = p — 1 + /, so (13) becomes 

(14) H , (z, </>) = \(P 2 Z + Pi) + \z 2 + h -4> 2 + l -c{z - f)4> -fz- f(l - /). 

where of course the constant term can be neglected. The following proposition 
makes explicit the integrals of motion in this case. 

Proposition 1. There exists a symplectic change of coordinates that trans- 
forms the Hamiltonian (14) to the form 

(15) H 0fi = Pi + Pi + + u\x\ 

PROOF. Perform the linear symplectic change of variables with parameter a 

z = x\ cos a — xi sin a, 

<f> = x\ sin a + X2 cos a, 

P Xl = P z cos a + Ptf, sin a, 

P X2 = —P z sin a + P^ cos a. 

Then naturally P£ t + P% 2 = P% + -Pj for any a. Imposing the condition that the 
mixed terms x\X2 vanish we get the choice 



-l + b± y/(6-l) 2 
tan a = 



Then the Hamiltonian becomes 



#o,o = \ {Pi + Pi + ujfxj + u 2 2 xl) 



where 



(16) ujf = cos 2 a + b sin 2 a + c cos a sin a = | ^1 + b + \/(— 1 + b) 2 + 

(17) io\ = sin 2 a + 6 cos 2 a - c cos a sin a = | ^1 + b - (-1 + b) 2 + c 2 ^j 



□ 
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4.3. The normal variational equations. We next calculate the full varia- 
tional equations of (10) along of the solution of (12), which can be written as 



(18) 



5p' 
59' 



= SP P , 

SP 9 
P 2 

= 5P d , 



2P e Sp 



8P' p = d J^-3^6p-5p-U<P-f5esm( 
' p 6 p 4 2 

8P' e = -fSp sin 6- fp56 cos 6, 



= -^p-bS<j>. 

It follows immediately from equations (18) that the tangencial equations along 
the invariant plane 6 = Pg = are 



Sp' = SP, 



SP 



4>i 



SP' = 



SPL 



Then this can be rewritten as 



/ Sp' ) 

5(f) 


( 


K 






\ 



--5p-b5(f>. 



1 



-§0 
-6 / 





( 5 p \ 

S(f> 




SP P 


) 


\ 6P* J 



The 2x2 lower matrix has precisely the eigenvalues —uj\ 2 given in (16), (17). Thus 
the eigenvalues of the full matrix are ±iu>±, H0J2 whenever w 2 2 7^ 0. The case 
w i 2 > i s obtained whenever c 2 < 46, and the general solution of this system is 

p(t) = A sin ui\t + B cosui2t, 
(f>(t) = C sin LU\t + D cos ix>it 

with arbitrary constants A, B, C and D. In the same way, the case uj\ 2 < is 
obtained whenever c 2 > 46, and the general solution of this system is 

p(t) = Ai sinhwit + B2 coshu^i, 
4>(t) = C 2 sinhwii + D2 coshc^i 



with arbitrary constants A2, B2, C2 and Z?2- 
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In general, for u)\ 2 7^ Oj provided whenever c 2 7^ 46, the general solution of such 
system is 

<f>(t) = C 3 e Wlt + D 3 e U2t 

with arbitrary constants A3, B3, C3 and D3. 

On the invariant plane the normal variational equation becomes 

86' = \sP e , 
P 2 

5P' = -fp56. 

Then we may rewrite this systems as a second-order equation of the form 

(19) 59" + ^-56' - -59 = 0. 

P P 

5. Galoisian analysis of the variational equation 

Considering the solution p(t) = Asin^ii + Bcosui2t, the variational equation 

(19) becomes to 

(20) d 2 9(t) + ^o S (u 1 t) -2Bu J2 Mu 2 t) _ f 

v 1 * w Asin(wii) +Bcos(w 2 t) Asin(wii) + 5 cos(w 2 t) 

where <9-t = ' = -4 (see appendices). To avoid triviality we assume / 7^ 0, cji 7^ and 
o; 2 7^ 0. By means of the change of dependent variable r\ = (As'mujit + B cosui2t)9 
this equation is transformed in the reduced form: 

(21) flfr(t) = f /-^^;nM)-i?^cosM)\ 
v ; * /w V Asin(cji<) + J Bcos(w 2 <) / 

The equations (20) and (21) cannot be algebrized using Hamiltonian Algebrization 
(see Appendix C) due to in general there is not exists r € Q such that L02 = ru>i, 
see [1, 2] and references therein. To apply Morales- Ramis theory, and in particular 
our result (Theorem 2.) we need one particular solution for the variational equation 
in where its Galois group is not virtually abelian. Thus, we can assume A = and 
the equation (21) is reduced to 

On the other hand, if we assume B = 0, we obtain 

< 23 » *<M^k;r4'' (t) ' 

Now, we proceed to apply the Hamiltonian algebrization for (22) considering B 7^ 0. 
We can see that x = x(t) ~ cos(oj2t) is a Hamiltonian change of variable, where 
yfa = — W2VI — x 2 . Therefore, the algebrization of (22) is 

(24) d 2 x fj(x)= ^-u^fftr), B x = y/Eid x , \=t ff( x ( t ))=r 1 (t). 
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Thus, the equation (24) is equivalent to 

(25) B^{x) + -^d x rj{x) + *~f X rf(x) = 0. 

x z — 1 uj^x[x a — 1) 

The equation (25) has four singularities: 0, 1, — 1, oo which are of regular type. This 
means that this equation can be transformed into a Heun equation: 

d 2 x y{x) + (! + -*- + -J-) d x y(x) + rt X - q y(x) = 0, 
\x x—1 x — aj x(x — l)(x — a) 

where e + 7 + <5 — ji — /3 = 1. With the translation i 4 i - 1 wc obtain the Heun 
equation 

„o^, , ( 1 1 \ „ ^, , —uskx + A + 1 . . , 

^ 26 ) + (to + 2^4) ^ + ^(x-l)( 8 -2) ^ = 0. 

with parameters 

a = 2, e = 7 = i <5 = 0, /x = 1, /? = -1, g= — ^ - 1. 

Through the change of variable £ = \J x(x — 2)rj we obtain the reduced form of the 
equation (26) to apply the Kovacic algorithm: 

„~ „ 3uj2 2 x 3 -UX + 9lu2 2 )x 2 + (8X + 3lu 2 2 )x + 3uj2 2 
(27) d 2 x i = rf , r = 



r = V r l = 0"2 = V?-oo = 0, aj = C*2 



4w 2 V (a; - 1) (x - 2) 2 

For our purposes, in order to apply the Kovacic algorithm (see Appendix B and see 
also [6]), we can write r as follows: 

_ -3 9u| - 8A A 3 8A + 9uf 

( ' T ~ 16(x - 2) 2 + 16w|(x - 2) + u? 2 {x - 1) ~ 16a 5 ~ 16w|x ' 

Now we start the analysis of the equation (27) applying Kovacic's algorithm. We 
can see that T = {0, 1, 2, 00} and o(r ) = °(r 2 ) = °(roo) = 2 and o(n) = 1. By case 
1, step 1, we fall into the conditions (ci) for c = 1, (c 2 ) for c = 0, 2 and (002) for 
00. In this way we obtain 

3 

4' 

- 1 + - 1 + 3 1 

a o= a 2=^ aj = a x =l, = -, = --. 

By step 2 we obtain D = {0} and by step 3 we obtain that Pq = 1 does not satisfy 
the relation (30). This means, for differential Galois theory (see Appendix A), that 
the Galois group of the variational equation is not a subgroup of the Borel group 
(except for A = or u; 2 = 0, avoided from the start). We follow with case 2 of 
Kovacic's algorithm in where we fall into the conditions (ci) for c = 1, (c 2 ) for 
c = 0, 2 and (002) for 00. In this way, by step 1, we obtain 

E =£ 2 = {1,2,3}, £i = {4}, £00 = {-2,2,6}. 

Now, by step 2, we obtain D = {0} and by step 3 we obtain that Po = 1 does 
not satisfy the relation (31). This means, again for differential Galois theory, that 
the Galois group of the variational equation is not conjugated to a subgroup of the 
infinite Dihedral group. Finally we look into the case 3 of Kovacic's algorithm. We 
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first consider m = 4, thus, by step 1, we fall into the conditions (ci) for c = 1, (02) 
for c = 0, 2 and (00) for 00. In this way we obtain 

E = E 2 = {3, 6, 9}, -Bi = {12}, ^ = {±6,0,12,18}. 

By step 2 we obtain D = {0} and for instance the rational function 9 and the 
polynomial S are given by 

6a; 2 -12a; + 2 3 2 

g » + 3»» + 2x ' S = - 3 + 3 - 2 + 2 - 

By step 3, we have the monic polynomial of degree 0, that is P = 1, and the 
sequence of polynomials is obtained by means of the relation (32): 



Pa = -1, 

P 3 = -6a; 2 + 12x - 2, 
P 2 = -6a; 4 + • • • + 3, 
Pi = 24a: 6 + ■ • • + 3, 



' 8 + --f, 

P_l = 3024a; 10 + • + 144 non-identically zero! 



P = -216a; 1 



and for instance the Galois Group is not the tetrahedral group. Now we consider 
m = 6, thus, by step 1, we fall into the conditions (ci) for c = 1, (02) for c = 0, 2 
and (00) for 00. In this way we obtain 

P = P 2 = {3,4,5,6,7,8,9}, E 1 = {12}, ^ = {±2^6,10,14,18}. 

By step 2 we obtain D = {0} and for instance the rational function 6 and the 
polynomial S are given by 

9a; 2 -18a; + 3 _ , 2 n 

By step 3, we have the monic polynomial of degree 0, that is P = 1, and the 
sequence of polynomials is obtained by means of the relation (32): 

Pe = -1, 

P 5 = -9.E 2 + 18a; - 3, 
Pa = -38016a; 12 + 

P_\ = 798336a; 14 + • • • + 6156 non-identically zero! 

and for instance the Galois Group is not the octahedral group. Finally we consider 
m = 12, thus, by step 1, we fall into the conditions (ci) for c = 1, (ca) for c = 0, 2 
and (00) for 00. In this way we obtain 

E Q = E 2 = {3, 4, 5, 6, 7, 8, 9} , E 1 = {12}, = {0, ±2, ±4, ±6, 8, 10, 12, 14, 16, 18}. 

By step 2 we obtain D = {0} and for instance the rational function and the 
polynomial S are given by 

18a; 2 - 36a: + 6 „ 3 2 



a; 3 + 3a; 2 + 2a; 



S = x A + 3x z + 2x. 
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By step 3, we have the monic polynomial of degree 0, that is P = 1, and the 
sequence of polynomials is obtained by means of the relation (32): 

Pl2 = -1, 

Pn = -18a; 2 + 36a; -6, 

P = -11280372719616a; 24 + • • • - 19737 f 8095 , 

P_i = 473775654223872a; 26 H h 85760073216 non-identically zero! 

and for instance the Galois Group is not the icosahedral group. In conclusion, the 
Galois group of the equation (27) is the connected and unsolvablc group SL(2,C). 

In the same way we can consider the case when c 2 > 46, that is, p(t) = 
A2 sinh uj\t+B2 coshu^t. Owing to the change of variable 1 1— > it transforms trigono- 
metric functions into hyperbolic ones and after similar change of variables as con- 
sidered before, we arrive again to a non-integrable Heun equation. In general, for 
c 2 7^ 46 which corresponds to p(t) = A^e^ 1 * + B^e U2t we arrive to a non-integrable 
confluent Heun equation in where not all singularity is of regular type. We can use 
the fact that the identity connected component of the Galois group is preserved 
under algcbrization process (see [1, 2]), thus, in our case, the Galois group of the 
normal variational equation will be the same no matter the way in where we alge- 
brize the differential equation. 

In this paper we shown that all the singularities of the normal variational equa- 
tion (27) are of regular type (Heun equation), as well we proven that such differential 
equation has not Liouvillian solutions and also we obtained that 



^ = |((i + &) + VFf+?), 

and wf = ± ((1 + 6) -7(1^6)2+^), 

being 

Thus, assuming fuiui2 7^ for any BeC*, we have proven the following result. 

Theorem 5.1. Let the parameters of the Wilberforce- spring pendulum satisfy 
c 2 — 46 ^ 0, then the system is not integrable through meromorphic first integrals. 

Proof. Observe that w 2 is always positive. On the other hand wf = if and 
only if c 2 = 46. □ 

Appendix A. Differential Galois theory 

The Galois theory of differential equations, also called Differential Galois The- 
ory has been developed by Picard, Vessiot, Kolchin and currently by a lot of re- 
searchers. In particular, we focus in the Galois theory of linear differential equations, 
also known as Picard- Vessiot theory. Following [2] and also [1], we present here an 
algebraic model for functions and the corresponding Galois theory. 

Differential Fields. Let K be a commutative field of characteristic zero. 
A derivation of K is a map d x : K — > K satisfying d x (a + 6) = d x a + d x b and 
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d x (ab) = d x a -b + a- d x b for all a, b £ K. We then say that (A, d x ) (or just A, when 
there is no ambiguity) is a differential field with the derivation d x . 
We assume that K contains an element x such that d x (x) = 1. Let C denote the 
field of constants of A: 

C = {ce K\d x c = 0}. 

It is also of characteristic zero and will be assumed to be algebraically closed. The 
coefficient field for a differential equation is defined as the smallest differential field 
containing all the coefficients of the equation. 

Due to we will mostly analyze second order linear homogeneous differential 
equations, i.e equations of the form 

C := d 2 x y + ad x y + by = 0, a 7 b £ K. 

so the rest of the theory will be explained in this context. 

Picard-Vessiot Extension. Let L be a differential field containing A (a 
differential extension of A). We say that L is a Picard-Vessiot extension of A 
for C if there exist two linearly independent y\ , y<x £ L solutions of C such that 
L = K{yx,y2) (i.e L = K(yi,y2,d x yi,d x y 2 )) and L and K has the same field of 
constants C. 

In what follows, we choose a Picard-Vessiot extension and the term "solution 
of £" will mean "solution of C in V . So any solution of £ is a linear combination 
(over C) of yi and yi- 

Differential Galois Groups A A'-automorphism a of the Picard-Vessiot ex- 
tension L is called a differential automorphism if it leaves K fixed and commutes 
with the derivation. This means that a(d x a) — d x (a(a)) for all a £ L and Va £ K, 
a{a) = a. 

The group of all differential automorphisms of L over A is called the differential 
Galois group of L over K and is denoted by DGal(A/A). 

Given a £ DGal(L/A), we see that {ayi,ay2} are also solutions of C. Hence 
there exists a matrix 

4,= (° J)eGL(2,C), 

such that 

o-((yi 2/2)) = (cr(yi) 17(2/2)) = (yi V2)A a . 

As a commutes with the derivation, this extends naturally to an action on a 
fundamental solution matrix of the companion first order system associated with 
C. 

- ( f/ 1 / 2 )) = ( m Vl \ m m \) = (V 1 / 2 ) A - 

\\o x yi d x y 2 J J \a(d x yi) a(d x y 2 ) J d x y 2 J 

This defines a faithful representation DGal(L/A') — > GL(2, C) and it is possible 
to consider DGal(A/A) as a subgroup of GL(2, C). It depends on the choice of the 
fundamental system {yi,y 2 }, but only up to conjugacy. 
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Recall that an algebraic group G is an algebraic manifold endowed with a group 
structure. Let GL(n, C) denote, as usual, the set of invertible n x n matrices with 
entries in C (and SL(n,C) be the set of matrices with determinant equal to 1). A 
linear algebraic group will be a subgroup of GL(n, C) equipped with a structure of 
algebraic group. One of the fundamental results of the Picard-Vessiot theory is the 
following theorem. 

The differential Galois group DGal(L/K) is an algebraic subgroup o/GL(2,C). 
In fact, the differential Galois group measures the algebraic relations between the 
solutions (and their derivatives) . It is sometimes viewed as the object which should 
tell "what algebra sees of the dynamics of the solutions" . 

In an algebraic group G, the largest connected algebraic subgroup of G contain- 
ing the identity, noted G°, is a normal subgroup of finite index. It is often called 
the connected component of the identity. If G = G then G is a connected group. 
When G° satisfies some property, we say that G virtually satisfies this property. For 
example, virtually solvability of G means solvability of G° and virtual abelianity of 
G means abelianity of G .. 

Lie-Kolchin Theorem. Let G C GL(2,C) be a virtually solvable group. Then 
G° is triangularizable, i.e it is conjugate to a subgroup of upper triangular matrices. 

Algebraic Subgroups of SL(2, C). We present below some examples of sub- 
groups of SL(2,C). 

Reducible subgroups These are the groups which leave a non-trivial subspace 
of V invariant. They are classified in two categories. 

Diagonal groups: the identity group: {e} = j ( q ^ ) \ , the n— roots: 



c 
c" 1 



1 



the multiplicative group: 



Triangular groups: the additive group: G a = ^ > ^ ^ ^ c?1 ^ uas ^~ 

roots: G*"} = {Yq c -ij ' c " = 1 ' d e C }' the Borel S rou P : 



( ° ) , 



I \ V o / 

Irreducible subgroups The infinite dihedral group (also called meta-abelian 
group): 

Doo = { f c^ 1 ) ' c£( C*|u|^ q J , d G C* | and its finite subgroups 

(where c and d spans the n-th roots of unity). 

There are also three other finite irreducible (primitive) groups: the tetrahcdral 
group A 4 L2 of order 24, the octahedral group Sf h2 of order 48, and the icosahcdral 
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group A5 2 of order 120. 

Integrability. We say that the linear differential equation C is (Liouville) 
integrable if the Picard-Vessiot extension L D K is obtained as a tower of differential 
fields K = Lq C Li C • • • C L m = L such that Li = Lj_x(r?) for i = 1, . . . ,m, where 
either 

(1) T) is algebraic over Li-\, that is ?y satisfies a polynomial equation with 
coefficients in L,_i. 

(2) ?7 is primitive over L,_i, that is c^Ty € Li— i. 

(3) ?y is exponential over L,_i, that is d x r\jr\ G L^_i. 

We remark that the usual terminology in differential algebra for integrable equa- 
tions is that the corresponding Picard-Vessiot extensions are called Liouvillian. The 
following theorem is due to Kolchin. 

TTie equation C is integrable if and only if DGb1(L / K) is virtually solvable. 



Appendix B. Kovacic algorithm 

Kovacic in 1986 (see [6]) introduced an algorithm to solve the differential equa- 
tion d^C = r(, where r G C(x). 

Each case in Kovacic's algorithm is related with each one of the algebraic sub- 
groups of SL(2, C) and the associated Riccatti equation 

d x v = t — v 2 = (y/r — v) (yfr + v), v = 

There are four cases in Kovacic's algorithm. Only for cases 1, 2 and 3 we can 
solve the differential equation, but for the case 4 the differential equation is not 
integrable. It is possible that Kovacic's algorithm can provide us only one solution 
(Ci), so that we can obtain the second solution (£2) through 

(29) C2 = Cl /f- 

For the differential equation given by 

d%( = r(, r = j, i,t6C[i], 

we use the following notations. 

(1) Denote by V be the set of (finite) poles of r, V = {c e C : t(c) = 0}. 

(2) Denote by T = V U {00}. 

(3) By the order of r at c € V , o(r c ), we mean the multiplicity of c as a pole 
of r. 

(4) By the order of r at 00, o (r m ) , we mean the order of 00 as a zero of r. 
That is o (7-00) = deg(rj) - deg(s). 

Case 1. In this case \\/r\ and [v^loo means the Laurent series of y/r at c and 
the Laurent series of y/r at 00 respectively. Furthermore, we define e{jp) as follows: 
if p G r, then e (p) G { + , — }. Finally, the complex numbers a+ , a~ , , a~ will be 
defined in the first step. If the differential equation has no poles it only can fall in 
this case. 
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Step 1. Search for each c G V and for oo the corresponding situation as follows: 
(co): If o (r c ) = 0, then 

[v^] c = 0, at = 0. 

(ci): If o (r c ) = 1, then 

[v^] c = 0, a± = l. 

(c 2 ): If o (r c ) = 2, and 

r = • • • + 6(a? — c)~ 2 + • • • , then 



[v^J c = 0, a c = . 

(c 3 ): If o (r c ) = 2u > 4, and 

r = (o (x - c)~" + ... + d (x - c)~ 2 ) 2 + b(x - c)- (l,+1 ) + • • • , then 

[Vr\ e = a(x - c)-' v + ... + d(x ~ c)~ 2 , a± = ^ (±^ + u ) ■ 

(001) : If o (t-qq) > 2, then 

[V^^^O, a+=0, a~=l. 

(002) : If ("Too) = 2, and r = • • • + bx 2 + • • • , then 

^00 = 0, a=o- ^ • 

(003) : If o (roo) = — 2« < 0, and 

r = {ax v + ... + df + bx"- 1 + • • • , then 

IVrL = ax " + ••• + rf ' and a °° = 2~V a _W 
Step £ Find D ^ defined by 

L> = |nGZ + :n = a<°°) - E a c (c) ,V( £ (p)) pgr | . 

If D = 0, then we should start with the case 2. Now, if Card(D) > 0, then for each 
n E D we search u> G C(x) such that 

w = e (00) [vf] ^ + J2 ( e (<o c + < {c) ( x - c ) _1 ) • 

cer' 

Step 5. For each nefl, search for a monic polynomial P n of degree n with 
(30) 5gP„ + 2ud x P n + (d x uj + uj 2 - r)P n = 0. 

If success is achieved then £1 = P n e$ ^ is a solution of the differential equation. 
Else, case 1 cannot hold. 

Case 2. Search for each c G V and for 00 the corresponding situation as follows: 

Step 1. Search for each c G V and 00 the sets E c ^ and Eoo ^ 0. For each 
c G r' and for 00 we define E c C Z and Eoo C Z as follows: 

(ci): If o (r c ) = 1, then £ c = {4} 
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(02) : If o (r c ) = 2, and r = ■ ■ ■ + b(x — c) 2 + • • • , then 

E c = 1 2 + kVTTIb : k = 0, ±2 j n Z. 

(c 3 ): If o (r c ) = u > 2, then £ c = {«} 

(c»i): Ifo(r 00 )> 2, then ^={0,2,4} 

(ooa)i If ( r oo) = 2, and r = • • • + bx 2 + ■ ■ ■ , then 

£oo = 1 2 + fcVl + 46 : fc = 0, ±2 1 n Z. 

(003) : If o (r^) = v < 2, then £00 = {v} 
Step 2. Find D ^ defined by 



^ feoo - ^e c V\/e p e s P! perl. 



If -D = 0, then we should start the case 3. Now, if Card(D) > 0, then for each 
n G D we search a rational function 6 defined by 

2 ^ X - c 

cer' 

Step 3. For each n £ D, search a monic polynomial P n of degree n, such that 

(31) d'iPn + WdlP n + (w x e + w 2 - 4r)d x p n + {die + wd x e + e A - Are - 2d x r) p n = o. 

If P n does not exist, then case 2 cannot hold. If such a polynomial is found, set 
4> — 9 + d x P n /P n and let w be a solution of 

uj 2 + cf)uj + - (d x <j> + 4> 2 - 2r) = 0. 

Then C,\ = e$ u is a solution of the differential equation. 

Case 3. Search for each cer' and for 00 the corresponding situation as follows: 

Step 1. Search for each c £ V and 00 the sets E c ^ and E^ ^ 0. For each 
ceT' and for 00 we define E c C Z and £"00 C Z as follows: 

(ci): If o (r c ) = 1, then E c = {12} 

(c 2 ): If o (r c ) = 2, and r = ••• + b(x — c)~ 2 + • • • , then 

12fc 



E c = {6+ VI + 46 : fc = 0,±l,...,± — , m = 4, 6, 12 ^ n Z. 

L to 2 

(00): If o (roo) = v > 2, and r = • • • + bx 2 + • • • , then 

Eoo = U+ — Vl + 46 : fc = 0,±l,...,± — , m = 4,6,12^n: 
m 2 



Step 2. Find L> 7^ defined by 



D = In e Z+ : n = — ^eoo - ^2 e cj ,Ve p £ E p , p £ T| . 

In this case we start with to = 4 to obtain the solution, afterwards to = 6 and 
finally m = 12. If D = 0, then the differential equation is not integrable because it 
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falls in the case 4. Now, if Card(D) > 0, then for each n <G D with its respective m, 
search a rational function 

~ 12 ^ x - c 

cer' 

and a polynomial S defined as 

S=U(x-c). 

cer' 

Step 3. Search for each n £ D, with its respective m, a monic polynomial 
P n = P of degree n, such that its coefficients can be determined recursively by 

, . P m = —P, and for i e {m,m — 1, 1,0}, 

[ ' Pi-i = Sd x Pi - ((m - *) d x S - SO) Pi-(m- i) (i + 1) S 2 rP l+1 , 

and for i = the polynomial P_i should be identically zero, i.e., P_i = 0. If P 
does not exist (P-i is not identically zero), then the differential equation is not 
intcgrablc because it falls in Case 4. Now, if P exists search oj such that 

i=0 v ' 

then a solution of the differential equation is given by 

C = e J , 

where w is solution of the previous polynomial of degree m. 



Appendix C. Hamiltonian algebrization 

In this section we follow [1, 2]. We recall that there are a lot of differential 
equations with coefficients that are not rational functions. For these differential 
equations it is useful, when is possible, to replace it by a new differential equation 
over the Riemann sphere P 1 (that is, with rational coefficients). To do this, we 
can use a change of variables. The equation over P 1 is called the algebraic form or 
algebrization of the original equation. 

Hamiltonian change of variable. A change of variable z = z{x) is called 
Hamiltonian if (z(x), d x z(x)) is a solution curve of the autonomous classical one 
degree of freedom Hamiltonian system 

d x z = d w H H = H(z,w) = ^ + V(z), 

o x w = —o z H v ' 2 

for some V € K, where K is a differential field. Thus, z = z(x) is a Hamiltonian 
change of variable if there exists a £ K such that (d x z) 2 = a(z). More specifically, 
if z = z(x) is a Hamiltonian change of variable, we can write d x z = yfa, which leads 
us to the following notation: d z = ^/ad z . 

We can see that d z is a derivation because satisfy d z (f + g) = d z f + d z g and 
the Leibnitz rules 



9z(f ■ g) = d z f ■ g + f ■ d z g, d z 



A d z f ■ g - f ■ d z g 



gj g 2 
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We can notice that the chain rule is given by d z (fog) = d g fogd z (g) ^ d g fogd z (g). 
The iteration of d z is given by 

3° z = i, d z = V^o z , = V^oJr 1 = V^d z (. . . (V^d z )) . 

•> v ' 

n times \/ad z 

We call Hamiltonian Algebrization to the algebrization process obtained by a Hamil- 
tonian change of variable. 

Hamiltonian Algebrization Theorem [1, 2]. Consider the systems of linear 
differential equations [A] and [A] given respectively by 

d x Y = -AY, d z Y = —AY, A=[a iJ ], A=$ ij ], Y=[ yil ], Y = [y a ], 

where G K = C(z(x), d x (z(x))), a,ij g C(z) C K = C(z,y/a), 1 < i < n, 
1 < J < ajj(x) = a,j(z(a;)) and yn(x) = yn(z(x)). Suppose that L and L are 
the Picard-Vessiot extensions of [A] and [A] respectively. If the transformation tp is 
given by 

X h-> Z 

then the following statements hold: 

(1) K~K, (K,d x )~(K,8 z )- 

(2) DGal(L/A') ~ DGal(Z/A) C DGal(L/C(z)). 

(3) (DGal(L/A))° ~ (DGal(i/C(z))°. 

A natural example of Hamiltonian Algebrization, and for instance of the in- 
troduction of the new derivative d z , is the case of second order linear differential 
equations. Consider d%y + ad x y + by = 0, using ip we obtain d z y + ad z y + by = 0, 
which is equivalent to 

(33) ad 2 z y + + yW) d z y + by = 0, 

where y(x) = y(z(x)), a(z(x)) = a(x) and b(z(x)) = b(x). 

In general, for y(x) = y(z(x)), the equation F{d x y, . . . , y, x) — with coef- 
ficients given by a% k {x) is transformed in the equation F(d z y, . . . ,y, z) = with 
coefficients given by at k (z), where ca k (x) = a,i k (z(x)). In particular, for y/a £ C(z) 
and a ik G C(z), the equation F(d™y, . . . , y, z) = is the Hamiltonian algebrization 
of F(d x l y, ...,y,x) =0. Now, if each derivation d x has order even, with a(z) and 
a,i k {z) being rational functions, then the equation F(d x l y, . . . , y, x) = admits an 
Hamiltonian Algebrization for z = z(x). An example, that illustrate this, is given 
by the following linear differential equation: 

d 2 x n y + a„_i(a;)<9^" _2 j/ + . . . + a 2 (x)d x y + ai(x)d%y + a (x)y = 0. 

Hamiltonian Algebrization Algorithm Theorem [1, 2]. In general is very 
difficult to find a suitable Hamiltonian change of variable, for this reason Hamil- 
tonian Algebrization is a method or procedure (not an algorithm!). For specific 
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families of differential equations we can obtain algorithms to apply Hamiltonian 
Algebrization, for example any differential equation 

F(d2y,d:- 1 y,...,d x y,y;e x ^,...,e XkX ) = 0, ti= Pi \i, p,e2, 1 < i < k, 

admits Hamiltonian Algebrization if and only if 

^eQ*, 1 < i < k, l<j<k 
A? 

and by means of the Hamiltonian change of variable z = e^ x we obtain the equation 

F(d?y,dr 1 y,---Ay,y,z) = o 

with coefficients in <C(z). 
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